RPA vs. exact shell-model correlation energies 
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The random phase approximation (RPA) builds in correlations left out by mean-field theory. In 
full Ohui shell-model spaces we calculate the Hartree-Fock + RPA binding energy, and compare it 
to exact diagonalization. We find that in general HF+RPA gives a very good approximation to the 
"exact" ground state energy. In those cases where RPA is less satisfactory, however, there is no 
obvious correlation with properties of the HF state, such as deformation or overlap with the exact 
ground state wavefunction. 
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I. INTRODUCTION 

One objective of nuclear theory, not reached yet, is an 
accurate method to compute the binding energy of nu- 
clides globally-that is, from helium to the superheavies 
beyond uranium, and to the driplines relevant to nucle- 
osynthesis. The most successful and widely referenced 
approaches are semiclassical, such as the finite-range 
droplet model jjj and the extended Thomas-Fermi plus 
Strutinski-integral model [Q. Exact, fully quantum me- 
chanical ab initio microscopic calculations starting from 
experimental NN scattering have been successful ^ [l| |f| 
but are computationally too demanding for all but the 
lightest nuclei. 

What about less ambitious but still microscopic ap- 
proaches? The method of choice for a global descrip- 
tion has been mean-field models such as Hartree-Fock 
(HF) with a phenomenological interaction (e.g. Skyrme) 
|6[. For computational simplicity Hartree-Fock is eas- 
iest to apply to spherical nuclei, although advances in 
computers allow one to treat deformed HF as well; in 
that case the deformed HF state breaks rotational sym- 
metry. Other symmetries that are broken by Hartree- 
Fock include translational invariance and (for Hartree- 
Fock-Bogoliubov) particle-number conservation. Mean- 
field theory thus leaves out both "kinematic" correla- 
tions, which restore broken symmetries and which must 
be corrected for approximately || , and "dynamical" cor- 
relations, which arise from the residual two-body inter- 
action. 

The interacting shell model (SM), which diagonalizcs 
an effective Hamiltonian in a basis of shell-model config- 
urations, includes correlations, respects rotational invari- 
ance and other symmetries, and yields accurate binding 
energies, but at a computational price. Shell model calcu- 
lations are currently limited to about 10 7 ~ 8 basis states, 
which restricts one to light to medium nuclei (A < 50) 
for full OTilo calculations. Furthermore, there is no global 
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shell model interaction; instead one must derive or fit an 
effective interaction for each regime such as the sd-shell 
(with valence nucleons between magic numbers of 8 and 
20) and the p/-shell (between magic numbers of 20 and 
40). 

The random phase approximation, or RPA, ||, ||] 
is an appealing middle ground because it adds correla- 
tions to the mean-field state and is numerically more 
tractable than full diagonalization in a SM basis; it also 
exactly separates out spurious modes induced by symme- 
try breaking. Because RPA implicitly treats 2-particle, 

2- hole correlations in the ground state, one might imag- 
ine RPA to be an approximation to a more sophisticated 
microscopic model, a kind of cleverly truncated interact- 
ing shell model so to speak. Unfortunately, this image 
is erroneous because, unlike HF and the interacting SM, 
RPA is not a variational theory, and in fact, RPA can 
overestimate the 'exact' binding energy. (This of course 
supposes one has the 'correct' effective interaction in the 
first place, a big question we will not attempt to answer). 

Surprisingly, tests of HF+RPA calculations of bind- 
ing energies in the literature have been scant and all 
have been against simple models: the Lipkin-Meshkov- 
Glick (LMG) model jl0[ Oj, a schematic interaction in 
a small shell-model space J12| ] (here HF+RPA did poorly 
in describing the "exact" binding energy), a generalized 

3- level LMG model |l3| and a two-level pairing model 
JLlj . Based particularly on the success of the latter two 
cases, Bertsch and Hagino advanced the idea of using 
the mean-field theory corrected through RPA for nuclear 
binding-energy systematics p5[ . 

This history motivates us to test the accuracy of 
RPA with respect to other microscopic calculations. We 
present a thorough study, comparing RPA to exact SM 
calculations in a large number of non-trivial systems. 

The paper is organized as follows. Section II presents 
the model space and the interaction used, shortly reviews 
the approach and compares the RPA and exact correla- 
tion energies for several nuclei. We find that RPA repro- 
duces well the ground state energy for a large number of 
nuclides. Because the HF+RPA occasionally misses the 
binding energy, we attempt in Section III to determine if 
the accuracy of the HF+RPA binding energy correlates 
with properties of the HF state, such as deformation ge- 
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ometry or the overlap of the HF state with the exact 
ground state. We find no such relationship. And be- 
cause RPA is often said to "restore" broken symmetries, 
in section IV we compare the kinetic energy of the RPA 
zero modes (associated with symmetry restoration) with 
the "kinematic" correlation energy obtained by project- 
ing the HF state onto good J. Finally, we summarize our 
conclusions in section V. 



II. CORRELATION ENERGY IN RPA 

In the interacting shell model one restricts the model 
space to a limited set of valence fermion states outside 
an inert core, interacting via a residual two-body Hamil- 
tonian. We consider nuclei throughout the sd-shell, with 
an 16 O core and using the Wildenthal USD interaction 
|L6| , and in the lower p/-shell, with a 40 Ca core and us- 
ing the modified KB3 interaction j[7). For each nucleus 
investigated in this paper we are able to perform a full 
Ohio diagonalization; in the same valence space we calcu- 
late the HF solution using separate Slater determinants 
for protons and neutrons. In this restricted model there 
is no mixing of radial degrees of freedom, the only de- 
grees of freedom being angular. We emphasize that this 
is nonetheless a much more sophisticated testing ground 
for HF+RPA than the usual toy models. 

Several textbooks discuss extensively the matrix for- 
mulation of RPA (see e.g., M), so we only summarize. 
One starts with a self-consistent mean-field solution, a 
Hartree-Fock Slater determinant, which can break sym- 
metries of the Hamiltonian and which ignores correla- 
tions. The correlation energy is the additional binding 
energy that arises from particle-hole correlations. For 
example, if one performs a shell-model diagonalization 
in the full space, then the shell-model correlation energy, 
which we consider to be the exact correlation energy, is 
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To derive the correlation energy in RPA, one treats 
the energy surface around the Hartree-Fock minimum 
as a multi-dimensional harmonic oscillator, which allows 
one to compute excited states as 1-particle, 1-hole (lp- 
lh) configurations, as well as to compute corrections to 
the ground state energy. This last is approximately the 
correlation energy from an admixture of 2p-2h configu- 
rations. Because in the mapping to a multi-dimensional 
harmonic oscillator one has bosonized the energy surface, 
RPA does not rigorously enforce the Pauli principle, so 
RPA implicitly assumes the admixture of 2p-2h states 
into the ground state is small. 

The resulting RPA matrix equation is well-known: 
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Most applications of RPA are to transition strengths 
and excitation energies, but one can also derive the RPA 
correlation energy: 



(4) 



One can also compute the correlation energy in the 
response-function formulation of RPA Jl^] . 

A particular feature of the RPA is the appearance of 
zero frequencies which correspond to conserved quanti- 
ties, sometimes interpreted as restoration of broken sym- 
metries H §]. In our restricted model space, the mean 
field can break the rotational symmetry; therefore, the 
number of zero RPA frequencies is related to the shape 
of the HF solution. Thus, when the HF solution is spheri- 
cal, there are no zero modes; when axially symmetric, we 
get two zero eigenfrequencies, while for triaxial shapes we 
obtain three. The zero modes require careful treatment, 
and one can rewrite Eqn. (0): 
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where Vp, are the canonical momentum operators corre- 
sponding to zero modes (for details, see [7j). Eqns. (Q) 
and (||) are equivalent, but the latter explicitly segregates 
the contribution from zero modes, that is, from restored 
symmetries @, |[ [l3[ |l4| . We compare the kinetic en- 
ergy due to the zero modes against the correlation energy 
obtained through projection of the HF state onto good 
J in section IV. 

Ideally one would hope for a good agreement between 
Erpa an d Es°m- Figures 1-5 compare the exact and 
RPA correlation energies for nuclides in the sd- and pf- 
shells respectively. Note that it is not our intention at 
this moment to test against the experimental data. In 
the sol-shell all the nuclei are numerically tractable for a 
full diagonalization; for an easier comparison we grouped 
into even-even (Fig. 1), even-odd (Fig. 2), odd-odd (Fig. 

3) and oxygen isotopes, that is, only neutrons active (Fig. 

4) . In the p/-shell the dimension of the SM basis rapidly 
increases with increasing the number of particles in the 
valence shell, limiting the number of exact cases we can 
hand; in Fig. 5(a) we present nuclides with both protons 
and neutrons in the active pf space, while Fig. 5(b) 
contains calcium isotopes (neutrons only). 

The general trend is a good agreement between E c ^[ 
and E R °p r A ; in each group, except for oxygen, there is 
only one or two nuclides where RPA significantly over- 
estimates or underestimates the correlation energy. In 
the so!-shell, if one excludes the oxygen results, the rms 
deviation for 41 nuclides is 870 keV; in the lower pf- 
shell, excluding calcium, the rms deviation is 480 keV 
for 11 nuclides. There was not a significant difference 
between even-even, odd-odd, and odd-even/even-odd nu- 
clides. The single-species results were significantly worse, 



FIG. 1: HF and HF+RPA energies relative to the exact SM 
ground-state energy, for even-even nuclides in the sd shell. 
The exact SM g.s. energy is set at zero, the upper bars are 
the HF energy and the lower bars the HF+RPA energy. 
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FIG. 2: Same as in Fig. 1 for odd-A nuclides in the sd shell. 



however: 1800 keV for 6 oxygen isotopes and 730 keV for 
7 calcium isotopes. We search for potential explanations 
for this difference, as described in section III below, but 
found none. 

Despite these limitations, based on the results pre- 
sented in Figs. 1-5, one can assert that RPA seems to be 
a promising candidate for a global microscopic theory of 
the ground state energies, as advanced in Jl5[. 



III. ANALYSIS OF RPA ACCURACY 

Despite our overall good agreement, RPA sometimes 
fails to accurately reproduce the exact correlation en- 
ergy. In this section we search for measures that corre- 
late with the success or failure of RPA, which may in 
turn point towards possible solutions. For example, if 
RPA fails to describe pairing correlations, one may turn 
to Hartree-Fock-Bogoliubov + quasi-particle RPA. For 
another example, one might expect RPA, which implic- 
itly assumes only small 2p-2h correlations, to be more 
accurate when the original HF state has a large overlap 



FIG. 3: Same as in Fig. 1 for odd-odd nuclides in the sd 
shell. 
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FIG. 4: Same as in Fig. 1 for oxygen isotopes (that is, 
neutrons only). 



with the exact ground state wavefunction. If this were 
the case, one might correct through some self-consistent 
or renormalizcd RPA scheme. 

To quantify our search use the ratio 



R = 



Tp cor r 
RPA 
rpcorr 



(0) 



R = 1 when HF+RPA agrees with the exact shell-model 
binding energy. 
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FIG. 5: Same as in Fig. 1 in the p/-shell for (a) nuclides 
with both protons and neutrons in the valence space and (b) 
calcium isotopes (neutrons only). 
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We begin by considering how well RPA works with dif- 
ferent geometries of the mean-field solution. In particu- 
lar, how well does the RPA correlation energy account for 
("restore") broken symmetries? Table I lists the defor- 
mation geometry (spherical, prolate, oblate, triaxial) for 
various nuclides; we compute /3, 7 by diagonalizing the 
mass quadrupole tensor of valence nucleons for the HF 
state |h| . R is close to one for a wide range of deforma- 
tion parameters. In fact, if anything RPA is less reliable 
for a spherical mean-field geometry, overestimating the 
correlation energy by as much as a factor of two. This 
result is surprising. In exact shell-model calculations, the 
excited states for these "spherical" nuclei are dominated 
by 1-particle 1-hole configurations, so one might expect 
RPA to be more successful than for "deformed" nuclei. 



TABLE I: Correlation energy ratio for select nuclides in sd 
and pf shells with different geometries of the HF solution. 



Nucleus 


P 


7 (degrees) 


E HF —E RPA 
E HF -Ef, M 


22q 







1.67 


24q 







1.83 


20 Ne 


0.46 





0.75 


22 Ne 


0.33 





0.97 


23p 


0.11 


60 


1.65 


24 Na 


0.24 


13 


1.17 


24 Mg 


0.29 


14 


0.92 


26 Al 


0.20 


33 


1.07 


28 Si 


0.20 


60 


0.90 


31p 


0.14 


40 


0.81 


32 g 


0.13 


32 


0.82 


34 g 


0.09 


48 


1.02 


34 C1 


0.10 


42 


0.91 


36 Ar 


0.09 


60 


0.81 


43 Sc 


0.28 


60 


0.80 


44 Ti 


0.44 





0.66 


45 Ti 


0.38 





0.97 


46 Ti 


0.35 





0.85 


46 Ca 


0.12 





1.04 


48 Ca 







1.66 



All of the spherical nuclei we found in the sd-shell were 
oxygen isotopes (valence neutrons only); this is not sur- 
prising as the proton-neutron interaction is well-known 
to induce deformation. To pursue this issue further, we 
took 28 Si, which normally has an oblate mean-field solu- 
tion, and lowered the (W5/2 single-particle energy, until 
the mean-field solution became spherical (filled d 5 / 2 for 
both protons and neutrons). In Fig. 6 we plot R as func- 
tion of the d§/2 single particle energy. Once again R is 
closer to 1 for the deformed regime than for the spherical, 
or closed-shell, regime. 

Now, as mentioned at the beginning of this section, 
one possible explanation for the poor results for spheri- 
cal nuclides is inadequate treatment of pairing. Because 
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FIG. 6: Phase transition in 28 Si as the CZ5/2 single-particle 
energy is lowered relative to the other single-particle energies. 
The solid line is R, the ratio of correlation energies, while the 
dashed line is the deformation /3. 



the onset of pairing is non-perturbative one might not 
expect RPA to describe pairing correlations well. In or- 
der to test this idea, we computed the expectation value 
of the pairing Hamiltonian in both the exact shell model 
calculation and in HF+RPA (2(J. We found, surprisingly, 
better agreement for oxygen and calcium nuclides and 
considerably poorer results for nuclides with more defor- 
mation and better RPA correlations energies. Hence we 
cannot point conclusively to pairing as the culprit. 

Another place RPA may stumble is in its implicit as- 
sumption that the admixture of 2p-2h correlations in the 
ground state is small. This shows up in two places, in 
the bosonization of the energy surface about the Hartree- 
Fock minimum, which does not enforce the Pauli princi- 
ple, and, in the computation of the RPA A and B ma- 
trices, the replacement of the RPA wavefunction by the 
HF state in equation (^) . To test this possibility, we cal- 
culate the overlap between the HF state and the exact 
wave function. This cannot be done by simply taking the 
projection 



O = (HF\SM). 



(J) 



In our SM calculations, we took the basis states to have 
J z = for A even and J z = 1/2 for A odd; the HF 
state is a superposition of states with good J z which are 
not restricted to the SM values. Therefore, any rotation 
of the Slater determinant can change the weight of J z 
components so that O is not invariant, although both 
Erpa an d Es°m are - We therefore define the overlap as 



O = Max({HF\SM}), 



(8) 



where the maximum is taken when rotating about the y 
axis. Note that when the exact ground state has J = 0, 
O is independent of orientation. 

Fig. 7 shows R as function of the overlap O; here we 
have included only nuclides in the se?-shell, but the results 
for pf are similar. O varies from close to zero to close 
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FIG. 7: The correlation energy ratio R vs. the overlap O, 
defined in Eq. (H), for select nuclei in sd shell. 



to remove contamination by spurious center-of-mass mo- 
tion. 



TABLE II: Kinetic energy of RPA zero modes compared with 
the correlation energy obtained after projection of the HF 
state on good J. 



Nucleus 


(MeV) 


KE (MeV) 


Eft? (MeV) 


20 Ne 


3.29 


2.55 


4.25 


22 Ne 


1.86 


2.29 


4.56 


28 Si 


2.67 


3.85 


6.97 


24 Mg 


5.34 


4.65 


7.10 


32 S 


2.35 


4.20 


8.08 


44 Ti 


1.18 


1.60 


3.72 


46 Ti 


1.82 


1.33 


3.24 



to one, but no compelling correlation with R appears. 
That is to say, the accuracy of the RPA correlation en- 
ergy does not seem predicated on how well the HF state 
approximates the exact ground state wavcfunction. 

A similar measure would be the amount of corre- 
lations in the RPA wavefunction, that is, N^f A = 
y\ l^mJ 2 - I n principle, RPA is more accurate when 
Npjf A is small, but, as in Fig. 7, we find no correlation 
between N^ A and R. Because of this, it is not not 
obvious that any renormalized, self-consistent, extended, 
etc., RPA scheme would in practice yield substantial im- 
provements. 



IV. KINETIC ENERGY OF SPURIOUS MODES 

One of the features that makes RPA appealing for cor- 
recting the mean-field results is its capacity to identify 
and separate out, as zero excitation energy, the modes 
associated with conserved quantities. Consequently, one 
can identify in Eq. (|^) a term in the correlation energy 
corresponding to the kinetic energy of the spurious solu- 
tion 



KEn = 



E 



2M„ 



(9) 



On the other hand, the simplest way to correct the bro- 
ken symmetries in the mean-field solution is to project it 
onto good quantum numbers. Table II presents the corre- 
lation energy obtained by projection of the HF Slater de- 
terminants onto good J, Ep°oj, and the kinetic energy of 
the spurious modes given by Eq. (Q); for comparison we 
have also included the exact correlation energy missing 
from the mean-field state. We conclude that KEq cannot 
be identified with the correlation energy obtained by di- 
rect restoration of symmetries obtained by projecting the 
HF state onto good quantum numbers. This may have 
relevance when if one goes to a multi-?iw space and wants 



V. SUMMARY AND OUTLOOK 

We have investigated the accuracy of HF+RPA for 
computing the nuclear binding energies, by comparing 
against the exact shell-model binding energies obtained 
through full OHlo diagonalization. In simpler models 

IT 
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|TJ], the RPA binding energy is generally, but 
not always fj^l , satisfactory. We found this also to be 
true with our significantly more complicated shell-model 
Hamiltonians. 

There were some non-trivial cases where RPA did not 
accurately estimate the exact binding energy; however, 
aside from the rather unhelpful observation that the in- 
accuracies tended to be worse for single-species cases, we 
could find no useful correlate with the accuracy of RPA. 
We did find, however, that RPA tended to yield very good 
results for deformed nuclides; despite this, however, the 
part of the RPA correlation energy identifiable with ro- 
tational kinetic energy had no relation with the energy 
gained from projection. 

How could one make RPA more reliable? There are 
two obvious steps. The first is to attempt to build pairing 
directly into the mean-field state through Hartree-Fock- 
Bogolyubov and then computing the quasiparticle-RPA 
correlations energy. Because we found no relation be- 
tween the accurate calculation of the expectation value 
of the pairing Hamiltonian and the accurate calculation 
of the total binding energy, however, we are far from con- 
fident that HFB+QRPA will yield significant gains. 

The other possibility is to correct RPA for its violation 
of the Pauli principle as well as using the RPA ground 
state wavefunction function in computing A and B ma- 
trices in equation (^). The literature documents numer- 
ous proposals along these lines: self-consistent renormal- 
ization of the mean-field single-particle energies and the 
residual interaction [|| |l| ; treating the particle-hole am- 
plitudes as variational parameters p2| ; rederivation of 
the RPA equations via the "number-operator method" 
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to account for ground state correlations |23|; and dress- 
ing the RPA phonons to account for ground state correla- 
tions 0] . None of these methods have become widely ac- 
cepted, and they overlap each other to a large degree, al- 
though the literature does not clarify similarities and dif- 
ferences. As in the case of ungarnished RPA, these pro- 
posals for improving RPA have only been tested against 



toys such as the Lipkin model Perhaps by imple- 

menting these various proposals for full-scale shell model 
calculations one could compare their efficacy and practi- 
cality. This we leave to future work. 
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